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ABSTRACT 


We report a sensitive search for redshifted CO (3-2) emission from the weak 
radio galaxy 53W002 at z=2.394. Maps at resolutions of 3" and 235 km s“^ show 
a signihcant emission peak within 0.5" of the optical and radio continuum peaks. 
The measured narrow band flux is approximately ten times the extrapolated 
cm-wavelength non-thermal radio continuum expected at 101.9 GHz and 
exhibits a spectral prohle implying a 540 km s“^ width (FWHM) at a systemic 
redshift z = 2.394 ±0.001 for CO (3-2). This emission has a total integrated flux 
of 1.51 ±0.2 Jy km s“^, approximately 4 times weaker than that previously seen 
in the lensed systems FSC10214±4724 and the Cloverleaf QSO. For a Galactic 
CO-to-H 2 conversion ratio, the implied molecular gas mass is 7.4x10^° Mq [Hq 
=75 km s“^ Mpc“^ and go =0.5). The CO emission is elongated at PA = 120° 
with a deconvolved major axis radius of 15 kpc (2.8"). This extension is along a 
similar direction to that seen in the cm-wave radio continuum and the optical 
but approximately three times larger. A velocity gradient is seen along the 
major axis, and if this structure is a (forming) disk, the implied dynamical mass 
is 9-22x10^° Mq at r < 15 kpc, assuming inclination i = 0° (edge-on). The 
magnitude of these masses and the similarity of the high gas-mass fraction are 
consistent with the host galaxy of 53W002 being a young galactic system but 
the metallicity (probably > O.IZq in order to produce the CO lines) implies 
signihcant heavy element production prior to z=2.4. This constitutes the hrst 
high redshift molecular gas which is detected in emission where there is probably 
no gravitational magnihcation. 

Subject headings: galaxies: individual (53W002) - galaxies: active - ISM: 
molecules — galaxies: formation 
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1. Introduction 


Early in the evolution of galaxies, there will be a phase in which the predominant 
mass component is interstellar and observations of this matter are critical to progress in 
understanding the initial building of galaxies and their stellar populations. Although the 
protogalactic gas should initially be atomic (or ionic), virtually all star formation models 
suggest that the gas becomes molecular during collapse to form stars. With relatively small 
metal enrichment (> 10% solar metallicity) in a prior, low-level starburst, this denser ISM 
becomes observable in the CO rotational transitions at mm-wavelengths. To date, there 
have been three conhrmed detections of CO emission at redshift z > 1: FSC10214-I-4724 
at z=2.28 ( Prown &: Vanden Bout 1992| , Solomon et ah 1992a| , [Scoville et ah 19^ ; the 
Cloverleaf QSO at z=2.56 (H1413-I-117 - Parvainis et al. 1994| , |Yun et ah 1997|) ; and 
BR1202-0725 at z=4.69 (Phta et al. 199^ , Pmont et al. 19"^) . In the hrst two (and 
possibly the last also, c.f. Pmont et al. 19961) , the high redshift system is gravitationally 
lensed by a foreground galaxy or cluster and the intrinsic (non-amplihed) properties are 
therefore too uncertain for broad inferences about protogalactic evolution. 


In this letter we report high resolution aperture synthesis observations designed to 
detect and image the molecular emission in the radio galaxy 53W002 and its associated 
galaxy group or cluster at Zopt = 2.390 (Windhorst et al. 1991, hereafter W91; Windhorst 
et al. 1994; Pascarelle et al. 1996a, 1996b). This cluster is clearly not subject to signihcant 
gravitational lensing based on the optical morphology and the fact that the separate galaxies 
have somewhat different redshifts and quite different spectra. 53W002 is unlikely to be 
lensed by the lower-luminosity z=0.581 elliptical galaxy that appears 5" to its North-West 
(c.f., W91; Windhorst et al. 1992, 1994). Moreover, no other lensed images at z=2.39 
are seen immediately surrounding this object in the narrow-band redshifted Lya WFPC2 
images of Pascarelle et al. (1996b). This radio-selected cluster was chosen for a sensitive 
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search for molecular emission since it appears to be a group of young, presumably gas-rich 
objects and Yamada et al. (1995) have reported a possible CO (1-0) emission feature here. 
Our observations clearly detect CO (3-2) emission (albeit at a level corresponding to a 
mass which is probably 5-14 times less than would be consistent with Yamada et al. ). 
Throughout this letter, we adopt Hq =75 km s“^ Mpc“^ and go =0.5. 


2. Observations 

53W002 was observed using the Owens Valley Millimeter Array in three conhgurations 
between February 1996 and January 1997. The array consists of six 10.4 m telescopes, 
and the longest baseline observed was 242 m. For 53W002 whose systemic optical redshift 
is z=2.390 (see W91), the CO (3-2) transition occurs at 101.95 GHz. The synthesized 
beams were 2.39x1.94" (uniform weighted) and 4.09x2.96" (natural), both at PA=—68°. 
The phase center, central redshift for the measured CO line, and adopted angular size and 
luminosity distances for 53W002 are given in Table 1. The primary beam of the 10.4 m 
telescopes is 70" and thus the observations cover all objects in the cluster within ~35", 
corresponding to 150 kpc at 53W002. Several other z=2.4 candidates are included in this 
area : 5, 11, etc. of Pascarelle et al. (1996a, b). The system temperatures were typically 
250-350 K in the signal side-band corrected for antenna and atmospheric losses and for 
atmospheric incoherence. Spectral resolution was provided by a digital correlator conhgured 
with 120x4 MHz channels (11.4 km s“^), covering a total velocity range of 1,300 km s“^. 
In addition, the continuum emission was measured in a 1 GHz bandwidth analog correlator. 
The total useable integration time on source was 61 hours distributed over 15 tracks in the 
three different telescope conhgurations. The nearby quasar 1828-1-487 (1.95 Jy at 102 GHz) 
was used to track the phase and gain variations, and Uranus (Tfe=125 K) and 3C 454.3 
were observed for absolute hux calibration. The positional accuracy of the resulting maps 
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is ~ 0.3". The data were calibrated using the standard Owens Valley array program mma 
( ^coville et ah 19^ and mapped using BITMAP ([Shepherd et ah 1994|) and the NRAO 
AIPS package. For the continuum maps (BW = 1 GHz), the rms noise was 0.20 and 0.48 
mJy for natural and uniform weighting, respectively. 


In order to predict more accurately the expected non-thermal radio continuum at 101.9 
GHz, we also obtained 25 min of data with the NRAO VLA at 14.9 GHz in the B-array 
on March 11, 1997. 3G 286 and 1658+476 were used for flux and gain calibration. The 
synthesized beam is 0.81 x 0.53 " (PA=~82° ) for these data. The source was unresolved 
(< 0.31" in diameter) with an integrated flux of 2.8 ±0.1 mJy. 


3. Results 


In Figure 1 the interferometric spectrnm integrated over the central 5" centered on 
53W002 is shown, smoothed to a velocity resolution of 235 km s“^ and sampled every 
half resolution element. The peak flux is 3.1 ± 0.5 mJy and the deconvolved line-width is 
540 ± 100 km s“^ (FWHM). For comparison, the peak flux of FSG10214±4724 is 14 ± 2 
mJy and the full linewidth is 250 km s~^( |l8coville et ah 1995| , Downes et ah l9^ ) and 
similar parameters hold for the GO emission from the Gloverleaf QSO ([Barvainis et ah 


1994|); thus the emission in the unlensed radio galaxy 53W002 is more than a factor of fonr 


weaker than that in the previously detected high redshift systems. 


In Figure 2 maps of the GO emission integrated over velocity ranges with AV=235 
km s“^ are shown. In Figure 3 the total GO line flux integrated over 200 MHz (590 km s“^) 
centered on the line is shown snperposed on the HST V-band image from Windhorst, Keel, 
& Pascarelle (1997). The flux in the integrated GO map is 1.51 ± 0.2 Jy km s“^ and the 
peak flux in the channel maps is 4.06 ± 0.89 mJy beam“^. No other GO emission featnres 
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were detected within our 70 " primary beam. We can therefore set upper limits at < 1.5 
mJy (2 (t) for the line flux at the positions of z=2.4 candidates 5, 11, etc of Pascarelle 
et al. (1996b). 

The deconvolved size of the emission feature seen in the integrated map (Figure 3) 
is 5.7 ± 1.3" by 1.7 ± 0.7" (FWHM) with major axis at PA=114° . At the angular size 
distance of 1.08 Gpc, the major axis diameter corresponds to 29.9 kpc, or a radius of 15 
kpc. In view of the limited signal-to-noise ratio in the maps, these size estimates are highly 
uncertain; however, it is clear that the emission is resolved since the centroid is shifted in 
the channel maps (see Figure 2) with positive velocity emission (relative to systemic) seen 
predominantly in the southeast and negative velocities in the northwest. The fact that 
this velocity gradient is along the major axis of the intensity distribution is suggestive of 
rotation, possibly in a forming disk. 

The peak flux in the narrow band maps (4 mJy) implies a beam averaged brightness 
temperature excess of 39 mK at A=2.9 mm or 0.13 K in the rest frame of the source. This 
corresponds to an absolute brightness temperature of 9.5 K in the rest frame after adding 
in the cosmic background radiation. No emission, above that expected due to the CO line, 
was detected in the 1 GHz (BW) continuum hlter at a level a = 200 /iJy - the CO line 
contribution averaged over 1 GHz is expected to be ~ 0.5 mJy. 

In Figure 4 the radio spectral energy distribution is shown, including the lower 
frequency fluxes of W91 and our measured flux at 15 GHz (2.8 ±0.2 mJy) and the measured 
continuum at 101.9 GHz corrected for the known line flux. A power law £t to the radio 
continuum (dashed line) gives a spectral index a = 1.2 ± 0.1. The line flux, averaged 
over 200 MHz (2.8 mJy) is shown as the solid dot to indicate that it far exceeds the flux 
expected from extension of the cm-wavelength radio continuum, and hence it must be 
line emission. We also note that there is no evidence at 15 GHz for a significant thermal 
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nuclear component setting in at higher frequencies since this point lies on the power-law 
extrapolation of the lower frequency data. Scaling Arp 220 for the H2 mass and luminosity 
distance, the expected thermal dust emission for 53W002 at 102 GHz (345 GHz rest frame) 
is about 0.1 mJy, and the free-free emission would be 10-100 times less (c.f. Scoville 
et al. 1997). 


4. Analysis and Discussion 

The center frequency of the GO(3-2) emission is 101.86 GHz, corresponding to a 
mean redshift of z=2.394±0.001. This is marginally different (350 km s“^) from the mean 
redshift z=2.390±0.001 for the optical/UV emission lines (W91, Pascarelle et al. 1996a, 
b). A similar offset with higher redshifts for the radio versus optical lines has been noted 
for lower redshift, high luminosity galaxies ( [Mirabel fc Sanders 1989]) . The offset might 
be explained if molecular gas and dust gas is falling into the galaxy, or if the optical/UV 
lines originate from gas flowing out of the galaxy with the symmetric redshifted emission 
extincted by dust within the system. 

From the integrated line flux of 1.51±0.2 Jy km s“^, we obtain a total GO (3-2) 
luminosity of L/iq=1.86x 10^° K km s“^ pc^, using equation (3) from Solomon et al. 
(1992b). (This is a factor of 3 weaker than that of FSG10214-I-4724.) For a Galactic 
conversion factor a=4 Mq (K km s“^ pc^)“^, the total molecular mass is 7.4x10^° Mq. 
The limits on the mass from other objects within the primary beam (eg. objects 5, 11, etc) 
are < 3x10^° Mq(2ct). If the clouds are hotter or denser than Galactic GMG’s, as would 
be expected in a vigorous starburst galaxy, the conversion factor should vary as jTco- 
To some extent, the decrease in the conversion factor due to the hotter gas may be offset 
by increased density. In Arp 220, an independent dynamical constraint on the conversion 
factor suggests a value ~ 0.45 times the Galactic value (jScoville, Yun, fc Bryant 1997] ). 






Based on the measured sizes and line-widths, it is also possible to estimate the 
dynamical mass associated with the CO emission feature. Assuming the galaxy is rotating, 
Mdyn = RVy 2 /Gsiv?i. For 1 / 1 / 2 , we adopt 200-250 km s“^, i.e. approximately half of 
the measured line width, and we take R= 10-15 kpc; we then obtain 0.9-2.2x 10^1 Mo 
for the galaxy, with no correction for inclination. The molecular gas mass, assuming the 
Galactic CO-to-H 2 conversion factor, would therefore constitute between 30-80% of the 
total dynamical mass. 


In present epoch galaxies, the gas mass fraction is typically < 10%; the probably much 
higher values indicated for 53W002 are consistent with this being a genuinely young galaxy, 
still in its initial phase of star formation. In nearby ultraluminous IRAS galaxies, high gas 
mass fractions are also indicated (cf. ^coville, Yun, fc Bryant 1997|) , but in these cases 
the gas is concentrated almost entirely in the central kpc - not over 30 kpc as in 53W002. 
Based on rest frame UV data for 53W002, W91 and Windhorst et al (1997) have derived 
star formation rate of 50-100 in order to account for the fluxes and colors. The age 

of the stellar population is estimated by them to be 0.3-0.5 Gyr in the center and 0.5-1 Gyr 
at ~ 10 kpc radius. At these rates, the estimated mass of molecular gas would be cycled 
into stars within ~ 10® yrs; this is also consistent with a young age for the system. 


53W002 also has a weak extension in the optical to the West 4 kpc; [Windhorst 


et al. 19^ which may be seen in the HST B- and V-band images (see Figure 3). The 


east-west orientation of this feature is similar to that of the radio source (cf. [Windhorst et 


al. 1991[) and the GO emission although the latter is extended on both sides of 53W002. 


Yamada et al. (1995) have reported detection of a possible GO (1-0) feature from 
53W002 at a central redshift z = 2.392 with peak and integrated fluxes of 5 mJy and 1.92 
Jy km s“^; however, it seems unlikely that this is the GO (1-0) counterpart of the GO 
(3-2) emission reported here in view of the much lower mass implied by our GO (3-2) flux 
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(and the marginally different redshift). For optical thick CO emission lines with the same 
brightness temperature in both CO transitions, the (3-2) flux should be 9 times greater (if 
the lines are optically thin, the ratio should be higher still). For low excitation GMC’s in 
the Galactic disk, the line temperature ratios are typically (1-0) : (2-1) : (3-2) = 1 : 0.7 : 0.4 
( panders et ah 1997|) . Thus even in the case of optically thick Galactic clouds, our line 
flux is approximately a factor of 5 weaker than what would be expected from the Yamada 
et al. feature, and in the more likely case of nearly equal brightness temperatures, it is a 
factor ~14 weaker. Lastly, it is possible that Yamada et al. have detected additional flux 
from the other 18 z=2.4 objects of Pascarelle et al. (1996b) in the vicinity of 53W002; 
however, our maps show no other comparably strong features within an area larger than 
their 50" diffraction beam. CO emission was not detected by us from the other objects 
in the interferometer field of view (70") at levels less than 1/3 of that seen in 53W002. 
These systems must have less massive ISM’s than the galaxy associated with 53W002. The 
CO luminosity of 53W002 is approximately three times that of the nearby ultraluminous 
IRAS galaxy Arp 220 (cf. pcoville, Yun, fc Bryant 1997 ) and is similar to the most gas-rich 
local systems such as IR14348-14 (panders et al. 19'MD and the radio galaxy PKS1345-I-12 
( [Mirabel fc Sanders 1989| , [Evans et al. 1996| ) . 


The large gas-mass found for 53W002 suggests that this is a very gas-rich and massive 
young galaxy - perhaps the progenitor of an elliptical or early-type spiral galaxy. In the 
course of its subsequent dynamical evolution, the gas in 53W002 may concentrate in the 
center (through dissipative energy loss as is likely to have happened in the nuclei of the 
nearby ultraluminous IRAS galaxies, cf. [Scoville et ah 1994[) . The central mass density 
could then become high enough to place the system in the area of the fundamental plane 
occupied by present-epoch elliptical galaxies ( [Kormendy fc Sanders 19^ . Already, the 
light profile in the envelope at radii 6-9 kpc apparently follows an law (Windhorst 
et al. 1992, 1994b, [Windhorst et al. 1997[) , for which W91 derived a luminous mass of 























2 — 4 X 10^^ Mq from nine-band photometry combined with its spectra. Lastly, we note 
that 53W002 is a compact steep spectrum radio source and at the present epoch, such radio 
sources are hosted usually by elliptical galaxies. 

In conclusion, we may be witnessing the birth of a luminous elliptical or early-type 
spiral galaxy. The copious amount of CO (and H 2 ) gas detected with OVRO suggest that 
such formation may have initially been associated with a massive disk, that is subsequently 
destroyed by mergers or tidal effects, eventually to leave a luminous bulge-dominated 
galaxy. 

The Owens Valley millimeter array is supported by NSF grants AST 93-14079 and 
AST 96-13717. The HST observations were supported by NASA grants GO-5308.*.93A and 
GO-5985.*.94A (to RAW and WCK) from STScI, which is operated by AURA, Inc., under 
NASA contract NAS5-26555. N. Scoville thanks Don Hall and Dave Sanders for hospitality 
at the Institute for Astronomy, University of Hawaii when part of this work was done. We 
also thank Aaron Evans for comments and suggestions on the manuscript. 



Table 1. Summary of 


RA (B1950) 

DEC (B1950) 

< 2 : >co 

Luminosity Distance^ 

Angular-Size Distance^ 

Mdynsin^i 

= 75 km s~^ Mpc“^, qo = 0.5 
'^using q;=4 Mq (K km pc^)“^ 
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Observations and Derived Properties 


17 ^ 12 ^ 59^86 

+50°18'51"3 
2.394 ±0.001 
12.41 Gpc 

1.08 Gpc (1" —5.24 kpc) 
(7.4 ± 1.5) X 10^° Mq 
(9 - 22) X 10^° Mq 
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Fig. 1.— Average spectrum of CO (3-2) emission in 53W002 smoothed to 80 MHz (235 
km s“^) resolution is shown for a 5" aperture centered on the position of 53W002. 

Fig. 2.— Maps of the CO emission in 53W002 averaged over 235 km s“^ velocity ranges are 
shown as a function of coordinate offset from the phase center. The cross marks the position 
of the 15 GHz continuum peak (see Table 1). The contours are —3, —2, +2, +3, +4, & +5 
times 0.89 mJy beam“^ (1<^)) and the natural weighted beam is shown. The velocities are 
relative to z=2.394. No signihcant continuum emission was detected (see §3). 

Fig. 3.— The CO (3-2) emission is shown integrated over all channels containing signihcant 
line emission (v=—270 to +270 km s “^). On the right panel, the 4 x 4" area immediately 
surrounding the radio continuum peak (cross) is shown zoomed in, with the CO emission 
contours superposed on the HST V-band image (Windhorst et al. 1997). The CO contours 
correspond to —3, —2, +2, +3, +4, and +5, times 0.44 mJy beam“^ (1<^)- 

Fig. 4.— The radio continuum spectrum of 53W002 is shown with points at 101.9 GHz to 
indicate the measured peak CO (3-2) line hux (solid dot = 2.8 mJy averaged over 200 MHz), 
the hux measured in the 1 GHz BW continuum hlter and the latter corrected for the GO 
line hux averaged over the 1 GHz hlter. The last, shown as the solid triangle, provides an 
upper limit to the true continuum at 101.9 GHz (i.e. 0.3+0.2 mJy). The radio huxes are 
from Windhorst et al. (1991) (600 MHz, 1.4 GHz, 8.4 GHz) and this article (15 GHz). The 
dashed line is a power law with a = 1.2 ht to the cm-wave data; the measured narrow band 
line hux at 101.9 GHz is far above the extrapolated radio continuum. 
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